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We report the existence of explicit flux configurations leading to stable de Sitter and 

Q H > Anti- de Sitter vacua, consisting on Type IIB compactifications on a 6-dimensional 

•^ anisotropic torus threaded with standard and S-dual invariant non-geometric fluxes 

p^j | in the presence of orientifold 3-planes. In all dS vacua the masses of the complex 

structure moduli are heavier than the Hubble scale suggesting that the axio-dilaton 

and Kahler moduli are natural candidates for small-field inflation. We also find that, 

at least for the reported solutions, compactifications on anisotropic torus increases 

M-r the number of moduli fields with sub-Hubble masses with respect to a compactifica- 

r- 1 ■ tion on an isotropic torus. 
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1 Introduction 

Over the last few years, the search for stable or metastable vacua has been the center of great interest by 
the string cosmology community [l|, |2j (for a recent review see [3j] and references therein). In particular, 
there is a class of models which has been studied in detailed concerning a Type II compactification on an 
isotropic or anisotropic torus in the presence of a Z2 X Z2 orientifold. The existence of stable dS vacua 
has been found by using different methods from algebraic formalisms to numerical calculations [4-I16J). 



In this work we continue our search, initiated in [17j |. for stable de Sitter (dS) and anti-de Sitter 
(AdS) vacua by implementing a genetic algorithm in a type IIB compactification on an anisotropic torus 
threaded with Ramond-Ramond (RR), Neveu-Schawrz-Neveu-Schwarz (NS-NS) and S-dual non-geometric 



fluxes denoted (Q, P) [18l . Il9l | in the presence of orientifold 3-planes. In this case we are not considering 



exotic orientifold planes or equivalently, odd- valued integer fluxes. 

Besides the necessity to incorporate the S-dual non- geometric fluxes P, we are interested in studying 
how their presence alters the existence of stable dS and/or AdS vacua and their influence on some cos- 
mological features as the existence of stable inflationary directions. We also are interested in comparing 
those same features in compactifications on isotropic or anisotropic torus, i.e., we want to study whether 
the torus anisotropy favors or not the presence of stable vacua and the existence of inflationary conditions. 

Although we are still far to present a full analysis, we found some generic features, as: 

1. There is always a hierarchy on the moduli masses and the complex structure is always heavier than 
the rest of the moduli . 

2. The supersymmetry (SUSY) breaking scale and the gravitino mass are very high, and their order of 
magnitude seems to be independent of whether the torus is isotropic or anisotropic as well as if we 
are considering both non-geometric fluxes Q, P or just one kind of them Q or P. 

3. It seems that the sub-Hubble moduli masses correspond to the Kahler moduli, favoring a possible 
inflation driven by them. In the isotropic case we find that only the Kahler moduli have sub-Hubble 
masses while in the anisotropic case, also the imaginary part of the axio-dilaton has a mass lighter 
than the Hubble scale. 

It is important to remark that our search is limited to some ad hoc selection of families of solutions 
in the flux configuration space which does not allow us to search the full space of flux configurations. 
Therefore, although the above features seems to be generic, they are the result of some concrete solutions 
we have found. A more complete analysis is still under work. 

Equally important is to notice that all these 4-dimensional models constructed from the inclusion of S- 
dual non-geometric fluxes are still not derived from a ten-dimensional string theory, although huge efforts 
have been made over the last few years. It is of great importance to see whether these models can be 
obtained from a compactification of a ten-dimensional theory 2(J. J21|. 



2 The effective scalar potential 



Following the analysis performed in [17| we look for stable de Sitter (dS) and Anti de Sitter (AdS) 
vacua in the context of Type IIB superstring compactification on an anisotropic torus threaded with 
non-geometric fluxes Q and their S-dual counterparts referred as P-fluxes. Standard orientifold three- 
planes acts by a Z2 symmetry on the internal space, i.e., we do not consider the presence of 07-planes 
or D-branes. The effective supersymmetric theory contains seven complex moduli fields, denoted by 



<&i = (pi + itp = {ti 2,3, S, C/1^,3}, where n is the complex structure, S the axio-dilaton and Ui the Kahlcr 
moduli. The corresponding superpotential is given by [l9( 

7 

W = P - $4^1 + Y, (^ " $4^3™) *m, (2-1) 

m=5 

where Pj = Pj(<I>i, $2, ^3) are polynomial of cubic order on the complex structure with real coefficients as 
follows (where the convention of summing on contracted indices is taken): 

fb(*l,2,3) = «00 - 4l*< + G(B*i*i " <$£$iQj$k, 
A ($1,2,3) = o? - 4i*i + o*2*i*i - o? 3 fc *i*j*fc, 
^($1,2,3) = a^ + a^^-a^^^-a^^^.^, 

Pf(*i,2, 3 ) = a$ + a$$i + <$*$& + <$»*&&, (2.2) 



where a™ 4 '' correspond to RR (r = 0), NS-NS (r = 1), Q (r=2) and P-fluxes (r = 3) integrated over 
3-cycles and with one leg over each torus, and the index m runs from 1 to 3. Notice that for the isotropic 
case the upper indices i,j, k are all equal since $1 = $2 = $3, while in the anisotropic torus, all of them 
are different. See Table [T] for our notation. 



Table 1: Integrated fluxes. Latin indices correspond to horizontal coordinates on the torus, while greek 
indices are related to vertical coordinates. Only invariant fluxes to the orientifold projection are shown. 
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Substituying the above polynomials in the superpotential W, we are able to compute the corresponding 
scalar potential given by 

v = e K (g i] g l g ] - 3) (2.3) 



where as usual Q %3 is the inverse of the Kahler metric, Q= = K — log(\W\ 2 ) and the Kahler potential K is 
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K = -J^H 2 ^)- (2.4) 



i=l 

The result is a scalar potential which is written as a function of the fluxes a™ 3 and 7 complex fields. 
In principle there is a huge number of flux configurations for which the scalar potential could have a 
stable minimum. However, the set of fluxes are constrained by tadpole conditions and Bianchi identities. 
Therefore, it is necessary to find a set of fluxes which fulfill these constraints before looking for extrema 
of the scalar potential. 

2.1 Constraints 

As it is known, there are basically two types of constraints a set of fluxes must fulfill. The first concerns 
the Tadpole conditions or roughly speaking the cancellation of D3-brane charges on T 6 : 



H 3 AF 3 = 16, (2.5) 

and the corresponding constraint on the Q and P -fluxes by imposing SL(2, Z) invariance and T-duality: 



Q-F 3 = I P-H 3 = (2.6) 

where as usual the Q and P fluxes contracts with a p-form to give a, p — 1 form. 

Similarly, the dual counterpart of Bianchi identity on H% given by dH 3 = is obtained by extending 
the algebra of isometry generators of the twisted background, where one can observe that the structure 
constants of the Lie algebra are in fact the non-geometric fluxes Q and P. The generalized identities are 



then given by (see |18l | for a formal derivation): 

Q[L Hmn\p - P[L F MN]P = 0, Q P Q R =0 (2.7) 

and 

Q [MN p L]P + p[MN Q L]P = Qj p [MN p L]P = Q (2 g) 

3 Stable dS and AdS 

By implementing a genetic algorithm we look for stable minima for the generated scalar potential by a set 
of fluxes satisfying the above constraints. Basically we concentrate our study to isotropic and anisotropic 
torus compactificationqfl. Such selection must be put by hand in our algorithm at the beginning of our 
calculations. See Appendix B for the algorithm code. 

3.1 AdS vacua 

We report the existence of at least 6 AdS vacua constructed from different flux configurations by compacti- 
fying on an isotropic torus. For the anisotropic case, we were not able to find a single stable case to report. 
Concerning the isotropic cases, we observe in general that the vacuum energy values are to high to be 
realistic. We show our results in Table [21 where we have omitted the indices related to the anisotropic case. 



isotropic compactifications with Q-fluxes in the absence of P-fiuxes were already studied in [lj 



Table 2: Flux configuration for AdS vacua. All cases correspond to compactification on an isotropic torus. 
The value for the scalar potential at its minimum is given in mass Planck units. 
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Notice that in the first three cases (A,B,C) there are not Q-fluxes, i.e., we are considering a compact- 
ification with only P-fluxes while the last three cases consider the presence of both non- geometric fluxes. 
The corresponding masses for all the 6 real stabilized moduli after diagonalization of the mass matrix (see 
[17l |) for details), are shown in Table [3l 



Table 3: Moduli squared masses for AdS vacua. 
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We observe that the energy vev and moduli masses seems to be independent of the type of fluxes 
present in the compactification. Also it is important to remark that the absence of anisotropic solutions 
might be a consequence of the method we are using here to find stable vacua. 



3.2 De-Sitter vacua 

We have found 4 explicit flux configurations fulfilling the tadpole and Bianchi constraints with a positive- 
valued minimum. We have not imposed extra conditions on the fluxes by the breaking of super symmetry, 
implying that supersymmetry is broken through all moduli. Similarly, in order to relax our system by 



avoiding more extra conditions, we have considered only even fluxes in our solutions, this means that we 
have not considered the presence of exotic orientifolds. 

There are two solutions which correspond to compactifications on isotropic torus with only P-fluxes 
(Case 1) and with both non-geometric fluxes (Q,P) (Case 2). We have also found 2 vacua related to 
compactifications on an anisotropic torus with only P-fluxes (Case 3) and with both fluxes (Q,P) (Case 
4) . Our results are shown in the Table UJ The value for the scalar potential at its minimum is given in 
mass Planck units. 



Table 4: Flux configuration for dS vacua. 
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Since we have not imposed extra conditions on the fluxes, we have assumed that SUSY is broken 
through all moduli, although the F-term related to the complex structure moduli has the higher value, 
indicating that SUSY is broken mainly by this complex field and therefore, establishing that the sgoldstino 
points towards this direction. For all cases, the F-term given by 



<F l >: 



=< e-^Thirje*' >, 



(3.1) 



are of order 10 3 in mass Planck units. Similarly, the gravitino mass m 3 /2 = e ' 2 W for the 4 cases is of 
order of 10 1 — 10 2 . On the other hand, the masses of all stabilized moduli are given in Tabled where the 



curvature contribution to the mass given by 



M 2 - 



D i d j r + -v x ij , 



(3.2) 



has been already considered as well as the required diagonalization of the mass matrix. In consequence, 
we can notice the following: 

• There is always a hierarchy on the moduli masses. Specifically the complex structure masses are 
always heavier than the rest of the moduli and they are upper-Hubble, implying that they cannot 
play the role of an inflaton in the small-field regime. The heaviest corresponding to the complex 
structure moduli, which breaks SUSY. Therefore we can notice that in all these models, the sgoldstino 
is almost orthogonal to possible inflationary directions (see below). 

• For the isotropic cases, there are subHubble masses corresponding to the Kahler moduli masses. 
This suggests the possibility that in case of having suitable conditions for small-field inflation, it 
would be driven by the Khaler moduli. 

• For the anisotropic cases we see that for both cases the sub-Hubble masses correspond to the Kahler 
moduli and to the imaginary part of the axio-dilaton. In this sense we can conclude that at least 
for these two cases, the chances to identify a possible slow-roll inflaton in the small-field regime 
increases. Wether this feature is generic in an anisotropic torus is still an open question. Notice that 
the value of the sub-Hubble masses of the Kahler fields come mainly from the curvature term in the 
mass formula 13.21 



Table 5: Moduli squared masses for dS vacua. 
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4 Final Comments 

In this work we report on the existence of AdS and dS vacua in the context of type IIB string compact- 
ification on isotropic and anisotropic torus threaded with RR, NS-NS and S-dual (Q, P) non-geometric 
fluxes in the presence of orientifold 3-planes. The case of P = on an isotropic torus was studied in [17|. 
By implementing a genetic algorithm we have been able to find 4 stable AdS vacua corresponding to a 



7 



compactification on isotropic 6-dimensional torus. For the anisotropic case, we were not capable to find a 
single stable AdS vacuum. 

Concerning the existence of dS vacua, we have found 4 different flux configurations leading to stable 
vacua. Although our algorithm does not allows us to search the entire flux configuration scape, we notice 
some generalities of the reported solutions: 

1. The complex structure masses are always larger than the axio-dilaton and Kahler moduli masses for 
both, the isotropic and anisotropic cases. The F-term related to the complex structure moduli is 
always larger than the rest, suggesting that the sgoldstino is directed along the complex structure 
derivative and since the masses are always larger than the Hubble scale, the complex structure is 
discarded as a possible inflaton in the small-field regime. 

2. The order of magnitude of the energy vev's, gravitino masses and F-terms seems to be independent 
of the isotropy or anisotropy of the 6-dimensional torus and of the presence of both non-geometric 
fluxes (Q,P) or just one kind of them Q or P. However we notice that the number moduli masses 
which are below the Hubble scale increases in the anisotropic case compared with the isotropic case. 
Therefore, at least for the cases we report, we observe that the size of sub-Hubble moduli increases in 
the anisotropic compactifications suggesting that in such scenarios the probability to find candidates 
for small-field inflation is also bigger, although deeper studies are needed. 

The models here presented corresponds only to a minuscule portion of the whole field configuration 
space. However we believe according to the form of the superpotential which clearly establishes a different 
role for the complex structure with respect to the other moduli, that the above features are rather generic. 
Clearly it is necessary to explore in more detail the space of solutions. It might be possible that corrections 
to the Kahler potential and/or incorporation of D-branes would change the values of masses and SUSY 
breaking scales to more realistic scenarios. 
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A Genetic code 

The following code shows the methodology employed to minimize the scalar potential for the anisotropic 
case. The remaining cases are taken as particular cases of this algorithm. The source code assumes 
that a scalar potential is stored on a handled function called V_anisotropic, whose entries correspond 
to the fluxes which satisfy both tadpole cancellation conditions and the Bianchi identities. The final 
points of interest given by the genetic algorithm are stored on a file called solucion_anisotropic.m and 
f lujo2_anisotropic.m. Once these variables are stored, the function "ToMathematica" translates all 
data into a Mathematica language. Finally, these solutions are further improved in mathematica by using 
the standard Newton method with an accuracy of 1000 decimals^ (the source code is not shown). 



clear all; 

clc 

\°/ up and lb contains the boundaries on the moduli space explored by the ga 

k=3; Y/oFamily solution 



We thank Saul Ramos- Sanchez for suggesting us this improvement in our method. 



lb=[.l, .0001,-500, .0001,-500,0.0001,-500,0.0001, . 
-500,0.001,-500,0.001,-500,0.001]; \°/ Lower bound 
Up=[lel0,lel0,lel0,lel0,lel0,lel0,lel0,lel0, . . . 
Iel0,lel0,lel0,lel0,lel0,lel0] ; \%Upper bound 
nvar=size(lb,2) ; \7oVariables number 
max=2000; \7oMaximum number of iterations 
time_ga_parallel=0 ; 
sol=zeros (max , 15) ; 
f luxes=zeros(max,64) ; 



fprintf('\n Minimization process begins\n'); 
for i=l:max 

[a0,all,al2,al3,a21,a22,a23,a3, . . . 

b0,bll,bl2,bl3,b21,b22,b23,b3, . . . 

C01,c02,c03,clll,cll2,cll3,cl21,cl22,cl23,cl31,cl32,cl33, . . . 

C211,c212,c213,c221,c222,c223,c231,c232,c233,c31,c32,c33, . . . 

d01,d02,d03,dlll,dll2,dll3,dl21,dl22,dl23,dl31,dl32,dl33, . . . 

d211,d212,d213,d221,d222,d223,d231,d232, . . . 

d233 , d31 , d32 , d33=constraint_solution(3 , 32) ; 
fluxes(i, :)=[a0,all,al2,al3,a21,a22,a23, . . . 

a3,b0,bll,bl2,bl3,b21,b22,b23,b3, . . . 

C01,c02,c03,clll,cll2,cll3,cl21,cl22,cl23,cl31,cl32,cl33, . . . 

C211,c212,c213,c221,c222,c223,c231,c232,c233,c31,c32,c33, . . . 

d01,d02,d03,dlll,dll2,dll3,dl21,dl22,dl23,dl31,dl32,dl33, . . . 

d211 , d212 , d213 , d221 , d222 , d223 , d231 , d232 , d233 , d31 , d32 , d33] ; 



°/ Here I call the function V_anisotropic.m which is a handle function 

V=V_anisotropicC(a0,all,al2,al3,a21,a22,a23, . . . 

a3,b0,bll,bl2,bl3,b21,b22,b23,b3, . . . 

C01,c02,c03,clll,cll2,cll3,cl21,cl22,cl23,cl31,cl32,cl33, . . . 

C211,c212,c213,c221,c222,c223,c231,c232,c233,c31,c32,c33, . . . 

d01,d02,d03,dlll,dll2,dll3,dl21,dl22,dl23,dl31,dl32,dl33, . . . 

d211 , d212 , d213 , d221 , d222 , d223 , d231 , d232 , d233 , d31 , d32 , d33) ; 
options = gaoptimset( 'Generations' ,500, 'TolCon' , le-6, 'TolFun' , le-10, 
'PopulationSize' ,200, 'MutationFcn' , Omutationadaptf easible, . . . 
'StallGenLimit' ,20, 'Display' , 'off , 'Vectorized' , 'on') ; 
startTime = tic; 

[x,fval,flag]=ga(V,nvar, [],[],[], [] ,1b, [] , [] .options) ; 
time_ga_parallel = toe (startTime) ; 
fprintf ('Iter . \°/ g takes \%g minutes. ' ,i,time_ga_parallel/60.) ; 



fprintf ('\n SolFluxes={a0->V/ d,all->\%d,al2->\7.d,al3->V/„d,a21->V/ d,a22->\7.d, 
a23->\7,d,a3->\7,d,bO->\7od,bll->\7od,bl2->\7d,bl3->\7,d,b21->\7,d,b22->\7.d, . . . 
b23->\7,d,b3->\7,d,c01->\7od,c02->\7od,c03->\7od,clll->\7od,cll2->\7.d,cll3->\7.d, . . 
cl21->\7.d,cl22->\7.d,cl23->\7od,cl31->\%d,cl32->\7,d,cl33->\7.d,c211->\7.d, . . . 



c212->V/„d,c213->\7.d,c221->\%d,c222->V/ d,c223->\7„d,c231->\7.d,c232->V/„d, . . . 
c233->V/.d,c31->V/ d,c32->\%d,c33->V/ d,d01->\7„d,d02->\ / d,d03->\7.d,dlll->V/ d,. . . 
dll2->\7.d,dll3->\7.d,dl21->\7od,dl22->\7od,dl23->\7,d,dl31->\7.d,dl32->\7.d, . . . 
dl33->\7.d,d211->\7.d,d212->\7od,d213->\7od,d221->\7od,d222->\7.d, . . . 

d223->\7.d,d231->\7od,d232->\7od,d233->\7od,d31->\7.d,d32->\7od,d33->\7od}; ' ,fluxes(i, :)) ; 
fprintf C\ntestsol={taurl->\7o0.25f ,tauil->\7»0.25f ,taur2->\7»0.25f ,taui2->\7o0.25f , . . . 
taur3->\7,0.25f ,taui3->\%0.25f ,Sr->\%0.25f ,Si->\7,0.25f ,Url->\7,0.25f ,Uil->\%0.25f , . . . 
Ur2->\7,0 . 25f , Ui2->\7.0 . 25f , Ur3->\%0 . 25f ,Ui3->\7,0 . 25f } ; \n ' , x) ; 
sol(i, : ) = [x,fval] ; 

end 

figure (1) ; 

loglog(l./sol(:,8),-sol(:,15),' x '); 

xlabel( 'g_s') ; 

ylabel('V_0'); 

title( 'Negative Lambda'); 

figure (2) ; 

loglog(l./sol(:,8),sol(:,15),' x '); 

xlabel( 'g_s') ; 

ylabel('V_0'); 

title( 'Positive Lambda'); 

7o Here I write the final solution ************************* 
save ( ' solucion_anisotropic ' , ' sol ' ) ; 
save ( ' f luj o2_anisotropic ' , 'fluxes ' ) ; 

ToMathematica 
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